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We calculate the dynamic single-particle and many-particle correlation functions at non-zero 
temperature in one-dimensional trapped repulsive Bose gases. The decay for increasing distance be- 
tween the points of these correlation functions is governed by a scaling exponent that has a universal 
expression in terms of observed quantities. This expression is valid in the weak-interaction Gross- 
Pitaevskii as well as in the strong-interaction Girardeau- Tonks limit, but the observed quantities 
involved depend on the interaction strength. The confining trap introduces a weak center-of-mass 
dependence in the scaling exponent. We also conjecture results for the density- density correlation 
function. 



I. INTRODUCTION 



Recent advances in experimental techniques have made it possible to study Bose gases in what are effectively two- 
dimensional and even one-dimensional (ID) traps [1], [2], [3], [4]. The confinement of the three-dimensional (3D) gas 
achieved by making the level spacing of the confining potential in one or two dimensions larger than the energy of 
the individual atoms. Because of the growing interest in coherent matter-wave interferometry and atom lasers, these 
recent developments have also revived theoretical interest in the properties of these effectively ID Bose systems. It 
is clearly important to understand in more detail the effects of confinement on, e.g., their ground-state and thermal 
properties including one-particle as well as many-particle correlations. It is well established [5] that a good theoretical 
framework for ultracold metal vapours is given by bosons with (in the case considered here) repulsive delta- function 
interactions. Wc can thus describe the ID system with a Hamiltonian 

7? = y jvi^ (x) {j^^ -I^ + V (x)^ ^ (x) + l^/it (x) (x) ^ {x) 4' (x)^ dx, (1) 

where ip{x) is a Bose field operator with its adjoint ip^x) and commutation relations [i^{x), ■4j^{x')] = 5{x — x'), etc., 
jjL is the chemical potential, g > is the interaction strength (the coupling constant), and V{x) = ^■VP'x^ is the 
harmonic trap potential. For V{x) = this system has been a prototype for exactly solvable models in ID [6], [7], 
[8], [9], [10], and many of its properties have already been worked out in great detail. 

In real physical systems the transition from 3D to ID behaviour has an interesting aspect in that the resulting 
effective density of atoms in ID can be either 'high' or 'low' depending on the parameters of the system. It is a 
feature of ID Bose systems that high density means weak interactions and low density strong interactions between 
the particles. This property opens up the possibility of realising experimentally the 'strong-coupling' limit [11] in 
which repulsive bosons display fermionic properties [6], [7]. This is the so-called Girardeau- Tonks regime. In this 
'true' ID regime in the sense of [12], the Bose-Einstein condensation in the sense of [12] is not realized, and the system 
is not coherent. 
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In contrast, in the limit of high effective ID densities of atoms, i.e., in the Gross-Pitaevskii, or the weak-interaction 
regime, short-range coherence builds up between the atoms, and the system forms at low enough temperatures a 
quasi condensate [13], [14], [15], [16], [17], [18], [19], [20]. True long-range order is in this regime destroyed by phase 
fluctuations, while density fluctuations are suppressed. It is thus evident that the local correlations between the 
atoms are difii'erent in these two regimes such that the 'local' two-point correlation function r(xi,a;i) tends to zero 
in the Girardeau-Tonks regime, and to unity in the Gross-Pitaevskii regime with quasi condensate [21]. At high 
temperatures or in the very small coupling-constant limit it obviously becomes two, the value for uncorrelated bosons. 
Of course we are not limited to two-point correlations: thus, e.g., the local three-point correlation function r(a;i ,xi,x{) 
is important because it is related to three-body recombination, an inelastic decay process [11]. Notice that in defining 
the full non-local correlation functions we need as many space variables being correlated because 

in the presence of a confining potential there is no longer translational invariance. Moreover at finite temperatures 
it becomes natural to correlate also time-like variables r for which the fields evolve in Matsubara representation 
'ip{x,T) = e~^^ip(x)e'^^ with H given by (1) and likewise for ■0'''(a;,r) [22]. Thus the two-point (single-particle) 
thermal correlation function which will be studied in this paper is given by 

r(a;i,ri;a;2,r2) = lTri)\xi,Ti)'ii>{x2,T2)^ , (2) 

and the n-point correlation function correspondingly. In these expressions is a 'time-ordering' operator which 
places the field operators from right to left in order of increasing r [22] . Note that these expressions like Eq. (2) are 
not translationally invariant in x but they are translationally invariant in r. Variables r lie in < r < /3, /3 = (fc^T)"-^ 
and T is temperature, and there is (for bosons) periodicity in r period /?. 

Despite the previous results for one or more space dimensions [15], [16], [23], [24], [25], [26], [27], [28] the decay 
of the single-particle and especially many-particle correlations in the presence of a confining potential have not been 
fully analysed. We address this problem for one space dimension in this paper, and point out that the way these 
correlations decay for increasing \xi — X2\ is governed by a universal scaling exponent that can be expressed in terms 
of observed quantities, and, at finite temperature, applies in both of the fundamental regimes described in [12]. This 
scaling exponent is analogous to the one that appears in the power-law decay of correlations in the homogeneous 
system of ID bosons at zero temperature. In the presence of a confining potential, this exponent becomes a function 
of the center-of-mass coordinate, but this dependence can be expected to be small. This result is valid in the weak- 
interaction Gross-Pitaevskii as well as the strong-interaction Girardeau-Tonks regime and thus straddles both of the 
two regimes distinguished in [12]. We also show how asymptotically all multi-particle correlations can be expressed 
in terms of the two-point correlations, and indicate how the results can be extended to density-density correlations. 
As all of these properties can be solved exactly by Bethe Ansatz techniques in the translationally invariant case of no 
confinement, this case provides a good point of reference, and we begin by reviewing it in the next section. 

Finally, the existence of the experiments like [1], [2], [3], [4], providing essentially one dimensional systems, means 
that it could become possible to check connected one dimensional theory in great detail at finite temperatures. A 
purpose of this paper is thus to point out that the methods employed here can be used to provide such a connected 
theory. 



II. SUMMARY OF EXACT RESULTS BY THE BETHE ANSATZ METHOD 

Let us consider the system described by the Hamiltonian (1). In the absence of the confining potential namely 
by setting V{x) = the Hamiltonian (1) becomes exactly solvable for its eigenstates and eigenvalues by the Bethe 
Ansatz method [9]. For a gas of N bosons on a ring of circumference L, namely for periodic boundary conditions of 
period L in the one dimensional space x, eigenfunctions of Hamiltonian (1) are [7], [9]: 



l$iv(Ai,...,Ajv)) 



7^/ 
/Mi 



dxi...dxNfN{xi,...,XN 1 Ai,...,Ajv)V'^(a;i)-V'^(a;jv) | 0), 



and these depend on A'' real parameters Ai, Ajv satisfying the so-called Bethe equations 



(3) 



=-n 



Aj ~\k + i{g/2) 
L^i A,- - A, - i{g/2) 



where g = {2m/h^)g is the renormalised coupling constant and 



l<k<j<N 



^_9^ if-{xj - Xk) 



2 (Ap - Ap J 



(4) 



(5) 
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Here the sum over all the permutations V of the numbers 1, 2, N is taken, e{x) = sign{x) is the sign function, and 
the factor is equal to 



n (A,- - Xk) 

j>k 



A^! n[(A, -Afe)2 + (5/2)2] 

j>k 



The Fock vacuum |0) is defined by the condition 'ijj{x)\0) =0. The corresponding eigenenergies are 



(6) 



The solutions Xj of the Bethc equations, Eq. (4), multiplied by h have a natural interpretation as the momenta of 
the particles each of mass m. The function /jv is a symmetric function of the variables xj and a continuous function 
of each xj. One can also see that /jv is an antisymmetric function of Xj. Hence, /jv = if Xj = Xk, j ^ k, in the 
whole coordinate space —oo<Xj < oo for each j. This property is the basis of the Pauli principle for one-dimensional 
interacting bosons. The wave functions (5) form a complete and orthonormal set in configuration space and thus the 
eigenfunctions (3) are orthogonal for the different sets of the solutions of the Bethe equations. At zero temperature 
(T = 0) the ground state of the system is defined by the solutions lying within the interval —Xp < Xj < Xp and these 
form a "Fermi sphere" in the momentum space. At an arbitrary value of the temperature (T ^ 0) the solutions are 
distributed along the real axis: — oo < Xj < oo. In the thermodynamic limit in which the number of particles N and 
the circumference of the ring L both go to infinity, N,L ^ oo, at a fixed density p = jj, the distribution function 
of the solutions of the Bethc equations is defined by the solutions of the Licb-Linigcr equations for T = [7] and by 
the solutions of the Yang- Yang equations for T > [8]. The free energy density £ of the gas is £ = limjv,L^oo En/L. 
The ground state with a particle density p > corresponds to positive values of the chemical potential in this exactly 
solvable case. 

In [9] it is proved that the correlation function (n < T2) 



{i!\xi,Tl)lp{X2,T2)} = 



_ tr (e-'3*V'^(xi,ri)^(x2,T2)) 



tr 



(7) 



may be expressed in the form 



{'tl>\xi,n)'tl){X2,T2)) 



{fir I ^'''(xi,Ti)V;(x2,r2) I I^t) 

{^T I ^t) 



(8) 



where | fix) is any of the eigenfunctions (3) contributing in the thermodynamic limit to the state of thermal equilib- 
rium. It depends on temperature through the distribution function of the solutions of the Bethc equations. For small 
nonzero temperatures and long distances \ xi — X2 \^ Ic (here Ic is a healing length Ic = h/^mpg) the asymptotics 
of the correlator (7) is 



{ll)\xi,Ti)lp{X2,T2)) 



\Xi - X2\+ ihv{Ti - T2)) 



1/9 ' 



(9) 



in which 6 is the scaling exponent 



2'K%p 

mv ' 



(10) 



and V is the sound velocity, = (p/m)^^. In this form the scaling exponent depends only on observable quantities, 
namely the density p, the sound velocity v, and the atomic mass m. The asymptotic result Eq. (9) with the scaling 
exponent (10) is valid for arbitrary values of the coupling constant < 5 < 00. For the so-called Girardeau- Tonks gas 
{g = 00) [6], [12], [29] Appendix B shows that 



6 = 2 



(11) 
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exactly as is well known [9]. For small coupling constants g we may express (10) in the form [31] 



= 27ThJ-^. (12) 
V 1T^9 

One may also express 9 (12) through the healing length l^: = 2'jTplc- In Appendix B we also give a simple derivation 
of these expressions (11), (12). 

Equation (9) means that the correlator (7) decays exponentially at long distances \xi — X2\ '> Ic > 

{il)\xi,Ti)ip{x2,T2)) ~ pexp ^-^^ ki - a;2|^ , (13) 



with the correlation length 



m 



(14) 



The result Eq. (9) also allows us to obtain the asymptotic behaviour of the correlation function at zero temperature. 
In the zero temperature limit (/3 ^ oo) the correlation function (9) transforms into the correlation function for vacuum 
fluctuations and has a power law decay 

{'4)'^{xi,Ti)i){X2,T2))T=0 ^ Y/9- (1^) 

\Xl - X2\+ ihv{Ti - T2) 

The distinction in the behaviour of the correlators at zero and finite temperatures may be explained by the 
formation of the Fermi sphere at zero temperature [9]. However this ground state of bosons cannot be considered as 
a true condensate since ti)V-'(.X2, ^2)) vanishes when \xi — a;2| — > 00 so that there is no long-range order in 

any general sense. As correlations decay algebraically, i.e., by a power law the system is scale invariant and can be 
thought to be at a critical point with 9 now a critical exponent. 



III. THE FUNCTIONAL INTEGRAL APPROACH 



A. The partition function 

We now evaluate strictly comparable results to those summarized in Section II by using methods of functional 
integration. We also obtain results which arc more general than those of Section II in that the effects of the confining 
potential V{x) = ^Q'^x^ are taken into account. 

The partition function Z for a one-dimensional gas with repulsive ^-function interactions of strength 5 > is given 
as a functional integral by [13], [30] 

Z= j e^I^'^lD^D^/;. (16) 
The classical action S of this system is equal to 

S'[V',^] = J dr J dx ^^jj{x,T)Kijj{x,T) — ^ijj{x,T)'ip{x,T)'ip{x,T)'ip{x,T)^ . (17) 

The differential operator in this action 5 is .ft' = — H, while the one particle Hamiltonian Ti includes the trap 

potential V{x) = and the chemical potential /i is also contained in H which is 7i = — 57^7^3- + V{x) — //,. 

We have introduced two c-number complex valued fields ip{x,T),'tp{x,T) and, consistent with the existence of the 
potential V{x), the boundary conditions for these fields are chosen to be vanishing at infinity for x (in the sense of 
quadratic integrability) and periodic, period /? = (kBT)^^, for r. These new c-mimber complex valued fields ip^ip are 
two independent fields, and they are introduced in a formal correspondence with the operator fields V^(x,r), tjj''{x,T) 
[13]. In (16), VipV-iij is the integration measure. 

For small enough temperatures we can expect to put ip{x, t) = i^oix, T)+ipe{x, t) and likewise for ip. The condensate 
variable ipoix^r) describes the condensate field and ipe{x,T) describes highly excited thermal particles. Also we can 
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expect to write the field variable tpo{x,T) = ipo{x) +^(x,t), where the field ipo{x) describes the ground state of the 
model at zero temperature (for the homogeneous case this field describes the particles in the Fermi sphere) , while the 
field ^(.T, r) describes the low lying excited particles (the phonon-likc excitations in the vicinity of the Fermi momenta 
in the homogeneous case). We shall also require that the fields tpo, fpe are orthogonal in the sense: / dxtpoix, r)V'e {x, r) 
= / dx'i()o{x,T)'tpe{x,T) = . Since by construction ifjo and tpe are orthogonal, and likewise for -tpo and 'tpe, we can 
write the integration measure as DipDip = Vipo'D'il^o'DiJe'D^e- 

We shall only consider terms in S up to quadratic (bilinear) in Ve,V'e- This means that we are making an approx- 
imation in which the non-condensed particles do not interact with each other. By doing this we can then actually 
integrate out the thermal fluctuations and obtain the effective action functional 



SefAi^o^o] = In y e^['^''^l2?V'e2?V^e, (18) 
depending only on the tpc'^Po variables. Then, 

Z= y e^«'"['^°''''"'lpVo2?V'o- (19) 
The effective action to lowest order in g is given by (see Eq. (A15) Appendix A): 

S'e//[V'o,'0o] = -PFncili) 

^ ( — S^ + K-V{x)\^bJx.T)-^^ 



+ J <^T-y dx|V'o(a;,r) ^ — ^ + A- y(a;)J ^/'o(a;,T) - -■(/;o(a;,r)V'o(a;,r)V'o(a;,r)'(/'o(a;,T)| , (20) 

where A = /i — 2gpnc{0) is a renormalized chemical potential and Fndjj) is the free energy of the non-ideal gas of 
thermal particles (see Appendix A). At this point it is reasonable to pass to new variables, namely the density p{x, r) 
and the phase ip{x,T) of the field iIjo{x,t). These variables are defined through, and compare [13], 



Mx, r) = ^^(^e'^^^'^\ Mx, t) = V^K^e-^^^^'"), 

where p and (p are two independent real fields. Now the integration measure 'D'tl)o{x,T)'Dtjjo{x,T) is changed to 
'Dp{x,T)T>(p{x,T). In these new variables the effective action (20) will be equal to 

Seff[p,v] = -0Fnc{lJ-)+i j dr y"fi.x|p9r<P+ ^9x(p5xV)| 

+ jyrjdx^^ {^pdl^p - pid^^f) + {A-V)p-^-p'Y (21) 

Here and below we denote the first order partial derivatives over r and x as dr and d^, respectively, and the second 
order ones as and d^. 

We consider the stationary phase approximation to Eq. (19). The appropriate extremum condition 6{Seff[p, <f]) = 
for the effective action (21) has the form of two thermal, r-dependent, Gross-Pitaevskii equations, namely 

idr^ + ^ (^^.'Vp - (d.^f) +{A-V)-gp = 0, 

-idrP+—d^{pd^ip) = 0. (22) 
The substitution of the solutions po, <^o of these equations into the effective action gives 

Seff[po,y:>o] = -l3Fnc{p) + ^ J dr J dxpl- (23) 

Here F„c (p) is the free energy of the non-ideal gas of thermal particles. The total free energy of the system is thus 

[16] F(M) = -r'^e//[po,H- 

At the Thomas- Fermi approximation [5], which is valid at low enough temperatures the kinetic energy term 
{d'^y/p)/ ^/p (so-called quantum pressure) in the first equation of Eqs. (22) may be omitted. The equilibrium solution 
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relative to the stationary ground state with drP = = dr^p is obtained by setting the velocity field v = m ^dx^p equal 
to zero, V = 0. The density profile is then given by, and compare with Eq. (A7) where ji is now replaced by A, 



p..(.).^p>.(.) = ^(i_^)e(i-g), (24) 



in which 6 is the Heaviside step function, and (recalling that V{x) = ^Q^x^) Rc is the radius of the ground state at 
zero temperature, R"^ = 

Following the original decomposition ipo{x,T) = ipo{x) + Ci^'-,''') wc suppose that the thermal fluctuations in the 
vicinity of the stationary state are small so that we may split the density profile as 

po{x, t) = ptf{x) + ttq{x, t). (25) 

The Gross-Pitaevskii equations (22) linearized around the equilibrium solution (po = Ptf,^ = const) will then take 
the form 

idrfo - QT^Q + -. dl-KQ = 0, 

4mpTF 



idr-Ko dx {pTFdxifo) = 0. (26) 

m 

Eliminating ipo and dropping the term proportional to h'^, we arrive at the thermal Stringari equation [32] which is 

^a%+a.((l-|)^=0, (27) 



where v is the sound velocity at the center of the trap. 



= = A. (28) 

m m 



The substitution tTq = e'^^'^u{x) transforms Eq. (27) into a Legendre-like equation 



U) 



Tx{{'-i)Tx<^^)=' 



^.^o«(2^) + -((l-^)-"(^)) =0- (29) 

Since the Thomas- Fermi profile (24) differs from zero only within \x\ < Rc, 'wc consider this equation for the fluctuating 
part u{x) also at \x\ < Rc- After the analytical continuation uj iE, the equation (29) possesses polynomial solutions, 
namely the Legendre polynomials Pn{x/Rc), if and only if 

Then this gives for the excitation spectrum En = hO.^ Mn±^ j33j 

B. The correlation functions 

In this Subsection we calculate the two-point correlation function for the non-homogeneous gas as defined by Eq. 
(2), 

r(a;i,Ti;a;2,T2) = (Trip\xi,Ti)tl>{x2,T2)) . (30) 



The case of multi-point correlators will be studied in Section V. 

We can express the correlator (30) as the ratio of two functional integrals [13], [34]: 

r(a;i,ri;a;2,r2) = Z-^ J e^^^'^H{xr,T^)i,{x2,T2)Vi,Vi,, (31) 

where the action S'[^, ^/i] is (17) and Z is the partition function Eq. (16). 
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We are interested in the long distance asyniptotics of the correlators when \ xi ~ X2 |^ Ic- The main contribution 
to the asymptotics is given by the low-lying excitations. By integrating out the fields tpg, 'ipf, included up to quadratic 
terms in Eq. (31), we find that the leading term of the asymptotics is equal to 

r(.„n;.„r2)^ , (32) 

where Seff[^oi ^o] is the effective action (20). We may rewrite this expression in terms of the density-phase variables 

/ exp \Seff\p,^] - iip{xi,Ti) +i(p{X2,T2) + ^lnp{xi,Ti) + ^lnp{x2,T2)]'Dp'Dlp 

r{xi,Ti;x2,T2) ^ f 7^ f iwr^Fi ' (33) 

J exp{Seff[p,^])VpVip 

where Sef/ip,^] is Eq. (21). 

At low enough temperatures we may change In p(a;i, n), In p(x2, T2) in (33) to \n pTF{xi),\ia ptf{x2) with ptf given 
by Eq. (24), since in the Thomas-Fermi regime density fluctuations are suppressed [35], [36]. The functional integrals 
in (33) can be evaluated by steepest descents [34], and the correlation function r(a;i,Ti;a;2,T2) can be expressed in 
the form 

r(a;i,Ti;a;2,T2) ~ exp^-S'e//[po, Vo] + Seff[pi,ipi\ 

- i(pi{xi,Ti) + iipi{X2,T2) + ^lnpTF{xi) + ^\npTF{X2)j, (34) 

where the term Seff[po, fo] is given by Eq. (23), and the fields po, (fo satisfy two stationary Gross-Pitaevskii equations 
(22). By definition, the fields pi,<^i are determined by the extremum condition: 



5 (^Seff[p,<f] - i(f{xi,Ti) + i(f{x2,T2) + ^lnpTF{xi) + ^ lnprF(a;2)^ = 0. 



(35) 



This variational equation leads to another pair of the Gross-Pitaevskii type equations. The first equation appears as 
a coefficient at the variation dp{x,T), 

idr'P + ^ - {d.V>f^ + (A - V{x)) -gp = 0, (36) 

while the second one is the coefficient at the variation 6ip{x,T), 

^2 

— idrp H dx ipdx^) = iS(x — xi )6(t — n) — iS(x — X2)6{t — T2). (37) 

m 

Note that ptf{x) is a fixed function, and it is not subjected to variation. Note also that the equation Eq. (37) is 
driven by (5-functions but Eq. (36) is not - corresponding to the suppression of density fluctuations. After substitution 
of the solutions pi, of Eqs. (36), (37) into the effective action Eq. (21), one obtains 

SeffiPl, ^1] = - 1 + f drj dxpl (38) 

We can furthermore consistently assume that, away from the boundaries, the density fluctuations are small so that 
pi{x, t) = pTF{x)+Tri{x, t), and tti is a slowly varying function. Consequently the terms y^d^y^ and dxiridxfi are 
small and can be dropped, and Eqs. (36), (37), when linearized around the Thomas- Fermi solution, can be expressed 
in the form 

Ti^ 2 

idripi ~ fl'TTi - — {dxVi) = 0, (39) 
^2 

-idrni-\ dxipTFdx'Pi) = i6{x - xi)6{t - n) - iS{x - X2)S{t - T2). (40) 

m 
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By differentiating Eq. (39) with respect to r, substituting the result into Eq. (40), and dropping the terms of order 
higher than g or , we find that 

1 

-d^ipi H {pTF{x)dxVi) = i5{x - xi)5{t - n) - iS{x - X2)6{t - T2). (41) 

g TTi 

It is convenient to rewrite this equation in the form 

-^—d^ipi + {pTF{x)dxipi) = i-2T{ d{x - xi)5{t - n) - ^(a; - X2)5{t - T2) \, (42) 

where v is the sound velocity at the center of the trap, Eq. (28), and ptf is defined in Eq. (24). The solution of this 

equation depends on the coordinates xi,ti,X2,T2 of the 5-sources so that (pi(x,T) = (/^i (.x, r; .xi, ti, X2, T2). Now, with 
the use of Eqs. (23), (38), (39), and (40), one can calculate the following contribution into the exponent in Eq. (34) 
which is 

-5e//[pO, fo] + Seff[Pl, Vl] 

- ^ ^'^ j ^^iPi~ Po) ^ ^ J J '^^ (/'I - Po) {Pi + Po) 

- 9 j dT JdxTnpo = J dr J dx ^idrifi - ^(dx(pij ^ po 

By substituting this expression into Eq. (34), we obtain for the correlation function 

T{xi,Ti;x2,T2) =i V Ptf{xi)ptf{x2) exp (^-^ (.^i(a;i,ri) - ipi{x2,T2))^, 

where ptf{x) is given by Eq. (24). 

It is convenient to express the solution of Eqs. (41), (42) in terms of the solution G{x,t;x' ,t') of the related 
equation 

^ dlG{x, r; x', r') + dJ (l - ^) dxG{x, r; x' ,t') \ = J-5{x - x')(5(r - r'). (43) 



(^)) 



and we call this equation as the '5-function driven Stringari equation'. This gives for the first order correlation function 
at xi ^ X2- 

r{xi,Ti;x2,T2) ~ V Ptf{xi)ptf{x2) (44) 

X exp (^-^ {G{Xi,Ti;X2,T2) + G{X2,T2;Xi,Ti)) + ^G{Xi,Ti;Xi,Ti) + ^G{x2,T2;X2,T2)y 

The function G{xi,ti;x2,T2) has the sense of the correlation function of phases, 

G{xi,Ti;X2,T2) = {ip{xi,Ti)(p{X2,T2)). (45) 

The substitution of this formula into (44) gives the well known result [13] 

T{xi,Ti;X2,T2) s/ PTF {xi)pTF {X2 ) exp (^-^{{0{xi,Ti) - 0{X2,T2)f)^. (46) 

Notice that the terms in the exponent of the formula (44) have different senses. The terms G'(a;i, ri; a;2,T2), 
G(x2,T2;a;i,Ti) in it depend on the relative position of the coordinates and thus define the long distance behaviour 
of the correlator. The terms G(a;i, ti; xi, n), G(x2, T2; X2, T2) each depend on a single set of coordinates and thus 
contribute to the amplitudes only. Thus we may express the correlation function in the form 

r(xi,Ti;X2,T2) ~ p{xi,Ti)p{x2, T2) exp {G(xi,Ti]X2,T2) + G{X2, T2] Xi, Ti))^, (47) 
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where ti), p(a;2, T2) are the renormahzed densities [16]. The sohition G{xi,ti:x2,T2) of Eq. (43) is defined up 
to an imaginary constant which has the sense of a global phase, and, as follows from (44), it docs not influence the 
phase fluctuations. 

The results of Section III are thus such that the excitation spectrum is given by the solution of the Stringari equation 
(27), while the correlation of phases {(p{x\,T\)(p{x2,T2)) is determined by G{x\,t\; X2,T2) which is the solution of the 
(5- function driven Stringari equation (43). 



IV. ASYMPTOTIC BEHAVIOURS OF THE CORRELATION FUNCTIONS 



A. Homogeneous gas 



Here we apply the method presented in the previous section, Section III, to the calculation of the asymptotics of 
the correlation function for the homogeneous gas. To describe the case when the trap is absent, V{x) = 0, we may 
equally send the radius of the condensate Rc to infinity. This way the Eq. (43) will take the form 

^ -d^^Gix,T;x',T') + dlGix,T;x',T') = J-dix-x')d{r-r'). (48) 



The solution of this equation is the correlation function of phases. We consider Eq. (48) in the ranges of arguments 
< a; < L and < r < /?. For | a; - x'\ < L, |t - r'| < /3, and /3"^ = ksT » {hv)/L, we obtain the solution 



TT 

sinh (\x — x'\ + ihv{T — r')) 



3 \---'\' (49) 



In the limit L ^ oo the second term in this expression disappears, and the substitution of the remaining expression 
in Eq. (44) gives the expected answer Eq. (9) for the asymptotics with the scaling exponent equal to 

^=^. (50) 



Using V = A/m and p = A/g for the sound velocity and density, respectively, one obtains the expression Eq. (10) 
for the scaling exponent. For small coupling constants we may re-express the result (50) in the form of Eq. (12). 

These asymptotics for the correlation functions are valid for the ring-shaped Bose-Einstein condensates, and the 
complete agreement between Eq. (9) found by the exact Bethe Ansatz method and the results of the functional 
integral methods in this case of homogeneous gas, confirms the validity of the functional integral method itself as 
described in its actual details above. And this should mean that the results achievable by the functional integral 
method in the inhomogeneous case V{x) ^ should be equally good. The case V{x) ^ cannot (so far) be treated 
by any Bethe Ansatz method as already mentioned. 

B. Trapped gas 

Let us consider the non- homogeneous Stringari equation (43) which determines through its solution G{x,t;x' ,t') 
the long distance behaviour of the correlation functions: 

^a^G(a;,r;a;^TO+a, (^(l- J)a,G(x,T;a;',r')) =^(5(a;-a;')(5(r-r'). (51) 

We are looking for the solutions which are finite, periodic in r, and quadratically integrable at \x\ < Rc- The 
substitution of the series 

G{x,t;x',t') = i^e-(---')G^(x,x'), (52) 
where co = n = 0, ±1, into Eq. (51) leads to an equation for the spectral density. 



10 



The solution of this equation is expressed in terms of the Legendre functions of the first and second kind, P^{x/Rc) 
and Qi,{x/ Re), which are hnearly independent solutions of the Legendre equation Eq. (29) with v taken to be 



We get: 



(54) 



where e{x — x') is the sign function: e(x) = sign{x). The validity of this statement can be checked by substitution of 
the result (54) into Eq. (53), and with the help of the equality [38] 

Q'Av)Pu{y) - QAy)K{y) = (i - vT'- 

For a stationary, r-independent correlator we have G{x;x') = ^Go{x,x'), where 

gRc 



Go{x,x') 



2% 

gR. 



^e(.-.'){Qo(^)~Qo(|;)} 
1 + ^) 



In 



(55) 



and we have used here the particular properties of the Legendre functions Pq{x) = 1 and Qo{x) = ^Inj^. This 
result leads to the following correlation function [15], [16], [27], [37] 



r(a;i;a;2) V Ptf {xi )ptf {x2 ) exp ( - In 

V 4/i 



\X1-X2\ 

2R^ 



(rr 1+3.-2 



1^:1-3:21 
2Rc 



(56) 



where is the scaling exponent as given in Eq. (10). 

To study the behaviour of the r-dependent correlation functions we must first notice that the Green's function of 
Eq. (53) is defined up to a solution of the homogeneous Legendre equation Eq. (29). To guarantee the convergence 
of the series (52), we may add such a term to the expression (54), and obtain thereby for the Green's function, with 

1^1 >0, 

^^(->-') = |^e(.-.')K(|;)P.(|;)-Q.(|^)P.(^)} 



gR. 



2h- 



(57) 



For \lo\ » hv/{2Rc), the asymptotics of the Legendre functions are [38]: 

2 



^-p; ) sin 

2i^sm^/ 



1^ + 



1. 



where cosO = x/Rc, < 6 < it, and | argz^j < ■k/2. Substituting these expressions into (57) we find the behaviour of 
the Green's function for large ji^j, 



G^{x,x') ~ - 



2hv\(j\ Vsineisin 



exp(-^\u;\\0-e'\). 



(58) 



In the quasi homogeneous limit, when |a;i — a;2| ^ ^1+^2 g^j^^j ^j^g coordinates xi and X2 are far away from the 
boundaries defined by Rc, so that jxi — X2\ <^ Rc, the leading term of the Green's function Eq. (55) is 
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Here S is the ccntcr-of-mass coordinate S = ^J^^, and v is the sound velocity Eq. (28). For the functions with 
oj ^ and at inverse temperatures larger than the lowest energy excitations, /? = {kBT)~^ ^ hv/{2Rc), we may use 
the asymptotics Eq. (58) to obtain 



A 



exp 



(-(M-' 



uj W X — X' 



■I) 



2fLVpTF{S) 

The substitution of Eqs. (59) and (60) in Eq. (52) gives now 

A 



G{x,t;x ,t ) 



■ln<^ 2 



inhvpTFiS) 

The correlation function for the Bose fields is in this case 



sinh — — - (la; — x'\ + iTivir — r')) 
ripv 



T{xi,Ti;X2,T2) 



\^Ptf{xi)ptf{x2) 



sinh ^ {\xi -X2\+ ihvin - T2)) 



i/e(s)' 



(60) 



(61) 



(62) 



where the scaling exponent 6{S) depends now on the center-of-mass coordinate S = ^^^^^ : 



0{S) 



2lTfipTF{S) 



mv 



Prom the result (62) it follows that the correlation function decreases exponentially at large distances, 

T{xi,Ti;X2,T2) :^ \/ Ptf{xi)ptf{x2) exp (- ^ '^^ 



2pn'^v'^PTF{S)J ' 



in which the correlation length depends on the center-of-mass coordinate and is 

^ 2h'(3prF{S) ^ hfJv 
m IT 

At zero temperature we find correspondingly that 

T{xi,Ti;X2,T2) = 



\/ptf{xi)ptf{x2) 



\X1 - X2\+ ihv{Ti - T2) 



(63) 



(64) 



(65) 



These expressions for the trapped gas are similar to those for the homogeneous gas, Eqs. (9) and (15), but the scaling 
exponents now depend on the centcr-of-mass coordinate. We thus find that our results for the trapped Bose gas 
precisely assume the forms of the exactly known results in the limit of vanishing trap potential. 

There is also another way to study the asymptotic behaviour of the r-dependent correlation function in the quasi- 
homogeneous limit [16]. We may suppose that the term {2x/B^)dxG{x,T;x' ,t') in Eq. (51) can be neglected, and 

that the factor [v^ jj-rpix)) in the first term in Eq. (51) can be substituted by the inverse squared sound velocity at 
the center of the trap. The corresponding equation for the correlation function of the phases in the quasi homogeneous 
limit is 



1 



:d^G{x,T;x',r') + dlG{x,T;x',T') 
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% h. V^pTF{^^) 

It is easy to check that the solution of this equation is Eq. (61). 



—J{x-x')5{t-t'). 



(66) 



V. MULTI-PARTICLE CORRELATION FUNCTIONS 

We consider finally the asymptotic behaviour of the multi-particle temperature-dependent correlation functions of 
the Bose fields i'^i', 



^m{xi,Ti; ...■,X2m,T2m) = {TT-lp\zi)...-lp\zrn)-lpi2.rn+l)---i^{z2m)), 



(67) 
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where z = (x,t). We show in particular that at low temperatures and for large relative distances Rij =\ Xi— xj \, the 
TO-particle correlation function can be expressed in terms of the two-point ones, 

^m{xi,n;...;X2m,r2m)- r{xi,Ti; X2,T2)~'"''\ (68) 

l<i<j<2m 

where li = 1 if i = 1, m, k = —1 if i = m + 1, 2m, and r(xi, Ti; Xj,Tj) = (TTip^{xi, Ti)tp{xj ,Tj)) is the two-point 
correlation function discussed in the previous sections. We suppose as well that the points Xi are lying inside the 
condensate and sufficiently far from its boundaries. The factorisation of the asymptotic behaviour of the correlation 
function Eq. (68) is similar to that in the non-trapped case [39] ; the lack of the translational invariance only introduces 
modifications to the two-point correlations. 

We can express Eq. (67) as the ratio of two functional integrals over the space of complex-valued functions (c.f. 
Section III): 

^m{Xl,Ti;...;X2m,T2m) = j (^^^*''^'^'4'{xl,Ti)..4{Xrn,rrn)i^{Xrn+l^'rm.+l)--4{^^ (69) 

where Z is the partition function Eq. (16). In order to estimate F^, Eq. (69), we proceed as in Section III. The 
leading term of the asymptotics is as follows: 

^ . ^ ^ /e^°//[^°'^°lV;o(zi)...V;o(z^)V>o(z^+i)...V;o(z2m)PV'oPVio 
r„(xi,ri;...;X2m,T2m)- J^sIJTvi^o ' ^ ' 

where Seff is the effective action Eq. (20). Repeating the arguments of Section III, we can rewrite the expression for 
the correlation function in the form 

(m 2m 2m \ 

Sef}\p,v]-iY. Vi7-k) + i E </'(Zfe) + 5 E ln/OTF(zfc)jI'pI'<p 

T. t \ k=l k=m+l fe=l /r7i\ 
^m{xuTl\...\X2m,T2m)'^ 7- 7^ ; ^ • (71) 

J exp [Seff[p,ip\)VpVip 

At low enough temperatures these remaining functional integrals can be evaluated by steepest descents, and the 
correlation function takes the form 

^m{xi,Ti; ...;X2m,T2m) — ^Xp^- Sgf f[Po, fo] + ^ef f[Pl, fl] 
m 2m ^ 2m \ 

-i^'Pi{^k) + iY^ (pi(zfe) -h -^lnpTF(zfc) j, (72) 

fc=l k=m+l k=l ' 

where the fields po) ¥'0 satisfy two stationary Gross-Pitaevskii Eqs. (22). By definition, the fields pi, <pi are determined 
by the extremum condition 

k—\ fc— m+1 h—\ 

and this variational equation also leads to a pair of the Gross-Pitaevskii type equations 

^2 2m 

- i drP H dx {pdx^p) =iy^ 5[x- Xk)S{T - rfe) - i 5{x - Xk)S{T - Tk). 

k—l k—m+1 

In the approximation considered in Subsection III.B (see Eqs. (39), (40)), these equations lead us to the following 
equation: 

^ m 2m 

-^—d^(p + dx {pTF{x)dxip) = i-^\Y\ S{x - Xk)S{T - r^) - S{x - Xk)S{T - Tk)\. (74) 

k—l k=m+l 
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One finds now that 

m 2m 
-Seff[PO,^o]+ Seff[pi,y:>l]^ ^(j^Vlizk) - <Pl(Zfc))- 

k=l /c=m+l 

Substitution of this expression in Eq. (72) leads to the correlation function 



2rn / . ra 1m >. 

r„(a;i,...,a;2m) ~ Jl VPTF(^exp (-^(^^iy9i(zfe) - ^ (^i(zfc)H, (75) 

k=\ ^ k=\ fe=m+l ' 



where the function i^i(z) = </5i(z; z^, "L-im) satisfies Eq. (74). We can then express the solution of Eq. (74) in terms 
of the function G(z, z'), which is the solution of Eq. (43), and this yields 



"K^'^^^^''^ ~ ^ </'i(zfe)) =-2!! H (G(zj,Zfc) + G(zfc,Zj)) 

fc=l fc=m+l j=l fc=m+l 

^ m— 1 m 2m — 1 lin ^ 2m 

+ 2^1 E (G(z,-,Zfe) + G(zfe,Zj)) + - ^ ^ (G(z,-,Zfc)+G(zfc,z,)) + -^G(z,-,z,-). (76) 

j=l k=j+l j=m+l *:=j+l i=l 

The first order correlation function at zi ^ Z2 is then equal to that of Eq. (44). Substitution of this expression in 
Eq. (76) gives the already indicated result Eq. (68): 

r„(zi,...,z2m) = n r~'-'^(Zi,Zj). 

l<i<j<2m 

VI. CONCLUSIONS 

It is plain that the functional integral methods used here to calculate correlation functions of ID bosons are very 
powerful. They extend the methods introduced by Popov [13], [34] to the case with a confining potential that breaks 
the translational symmetry of the system. A key new ingredient introduced here is the Green's function G that 
satisfies the (5-function driven Stringari equation Eq. (43). This Green's function was shown to be precisely the 
correlation of phases. We thus obtain the results quoted in [15], [16] in a very natural way, for the Green's function 
could then be used to obtain the two-point correlation function. Note how density fluctuations were suppressed by 
the step at which ^- function driving terms were deliberately neglected from the Eq. (39). It is rather remarkable 
that for a vanishing trap potential the two-point correlation function found by the approximated functional-integral 
method coincides precisely with the exact Bethe Ansatz result, especially when one takes into account the changed 
boundary conditions. 

The 'universal' feature in the two-point correlation functions obtained with or without a trap potential is the 
scaling exponent 6 that governs the decay of correlations for increasing separation of the two points. At non-zero 
temperatures 9 is proportional to the correlation length and at zero temperature it becomes a critical exponent in the 
algebraic decay of correlations. We express this scaling exponent in terms of the observed quantities, sound velocity 
V and the density p of the quasi condensate, which means we need not calculate the T-matrix as, e.g., in [24]. It also 
means that the expression we find for 6 is valid for all coupling strengths: the two physical quantities involved do 
depend on coupling strength so that = 2 in the exact Girardeau- Tonks limit g oo, and > 2 for any finite g. In 
the presence of a confining potential the decay of correlations is formally exactly the same, but the scaling exponent 
will now depend on the center-of-mass coordinate, 

^ 2TrnpTF{S) 

mv 

with S = ^J-i^. As we expect the S'-dependence of d to be small, it should approach two in the Girardeau- Tonks 
limit also in this case. In the weak-interaction Gross-Pitaevskii limit d can also be expressed in the form 

Q[S) = 27rn. 



mgpTF{0)' 
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As all quantities that enter these expressions for 9 are themselves observables, a measurement of the decay of cor- 
relations in a trapped Bose system would provide a sensitive test of the theoretical framework based on repulsive 
5-function interactions between the particles. With a sensitive enough measurement one should also be able to detect 
the effects of confinement, i.e., the dependence on the center-of-mass coordinate S in the decay of the correlations. 

As experiments on trapped Bose condensates have become increasingly more sophisticated, a need for better under- 
standing of also many-particle correlations has become evident also. So far little has been known about the 2m-point 
correlation functions for m > 1 . One of the important results of this paper is thus the demonstration that all 2m-point 
correlations in trapped Bose systems can be expressed asymptotically entirely in terms of the two-point correlations. 
Note that, based on the results given here, and the exact results wc have reported earlier on a related Bose system 
[31], we can also propose that the asymptotic form of the density-density correlator for trapped Bose systems as 

Tr1p\xi,Ti)'4>{xi,Ti)'lp'' {X2,T2)'tp{X2,T2)'j - Ptf{xi)Ptf{X2) 

~ ^ + ^) + COS [27r(a;i - X2)ptf{S)] 

with w = sinh{(7r|a;i — X2\/h(}v) + ihv{Ti — T2)}. Here A and B are real constants. At large separations of points 
this density-density correlator decays as Ae~'^'^^~^^^ with the correlation length = {2'k)~^%I3v. Notice that the 

scaling exponent 6 appears here in the oscillating term which is not the leading term asymptotically. One would nev- 
ertheless expect that the oscillating behaviour in the asymptotic density-density correlations should be experimentally 
observable. 

Finally we would like to emphasize that the methods applied here are not restricted to ID but can equally well be 
applied in 2D and 3D [16], [42]. These methods are not restricted to the calculation of correlation functions, and can 
be used to calculate all of the thermodynamic properties of interacting Bose systems, and can be extended beyond 
the Thomas-Fermi as well as the mean-field approximations. 
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Appendix A 



Hero we shall give a formal derivation of the effective action S'e//[f/'o, fAo], Eq. (20). We use the field-theoretical 
approach of loop expansion [40]. Following the general prescription of Section III, we split the initial fields V'l V' iiito 
ip = "^0 + V'e) f/' = '0O + V'e SO that the action Slip, ^j], Eq. (17), is divided into three terms: 



where Scond is the action of the condensed particles in the so-called tree-level approximation [43] 

r0 



The action of the highly excited thermal fluctuations of the non-condensed particles is given by 



SfreebPe,'4'e] = \ j ^'^ j 



dx{lpe,1pe)G~ 



and Sint is the part of the action that describes the interaction of condensed particles with the thermal ones. 



dr / dx 



[ipe{x,T) 



tpo{x,T) -I- V'e(a;,T) 



K- - gipoipo 



In Eqs. (A2)-(A4) we have defined the differential operators. 



n = 



2m dx^ 



M + y(x), 



(Al) 



S..U^.,n ^ ldrJdx{Mx,r)K,Mx,r) - f ^„(.,.)^„(.,.)^„(.,.)^„(.,.)} . (^2) 



(^3) 



Mx,t)}. {AA) 
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and the matrix operator 
where 



= Go ^ - S, (A5) 







The approximation used here for Sint, Eq. (A4), imphes that the non-condensed excitations do not interact with each 
other. 

For our purposes it is enough to consider the functional integral in leading order as given by the method of stationary 
phase (steepest descents). Therefore we choose 'ipci'o as the stationary points of Scond (A-2). They are governed by 
equations of the Gross-Pitaevskii type, 

d d"^ \ 

d^^^d^^^~ ^^^^ I ~ 9{^o^o)i^o = 0, 



As soon as ^/jq. '^^^ the solutions of these equations, Sint drops out of Eq. (Al), and the dynamics of ipe^ ipe is just 
given to lowest order by Sfree, Eq. (A3). The latter depends non-trivially on ■ipo, i^o through S in G~^, Eq. (A5). We 
denote a r-independent solution of the Gross-Pitaevskii Eqs. (A6), written in the Thomas-Fermi approximation as 

^oV'o = Ptf{x; m) = ^ (m - V{x)) e (/X - V{x)) , {A7) 

where one should notice the dependence of the Thomas-Fermi profile on the 'bare' chemical potential /x. The inte- 
gration over V'eiV'e in Eq. (18) is now Gaussian, and gives a one-loop corrected effective action [13] in terms of the 
condensate variables tpci^o only, 

Seff[i^o,i>o] = Scond[i^o,^o] - ilndet(G-^). (AS) 

Here is the matrix operator (A5). 

To assign a meaning to the resulting expression for the effective action (A8), we have to regularize the determinant 
det(G~^). Our is already written in the form of a 2 x 2 matrix Dyson equation (A5), where the entries of the 
matrix S('!/'o) ipo) play the role of the normal and anomalous self-energy parts. The Dyson equation defines the matrix 
Green function G = G{ijjo,i>o) {i-e., the propagator matrix) of the fields tpe,4'e- The matrix G is given by the formal 
inverse 



G=(Go-1-s) \ {A9) 



The matrix operator G^^ (A5) can be formally diagonalized by means of the famous {u, ?;)-transformation of Bo- 
goliubov [41], [35]. Compatibility of the corresponding equations for the unknown fields u,v defines a quasi-classical 
spectrum of the elementary excitations [32]. 

For our purposes it is appropriate to represent G~^ in the form 

G-i = Go^-± = g-^ - (S - 2gpTF{x; m)/) , (^10) 
where J is a 2 x 2 unit matrix, and is defined as 

Here ptf{x',h) is the solution (A7) of Eq. (AG), and the expression (AlO) implies that we have simply added and 
subtracted 2gpTF{x] p) in the diagonal of the matrix operator G~^. A formal inverse of can be found from the 
equation which defines the corresponding Green's functions Q±: 
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Applying the relation Indet = trln, which is valid for both matrices and operators, we obtain 

-ilndetG-i = -^trln (/ - g{t - 2gpTF{x; ^i)!)) - ^Indet ( ^+ ^_ \ . {A12) 



The first term on the r.h.s. is free of divergences and now we have to assign a meaning to the infinite-order determinant 
on the r.h.s. of (A12). The operators IC± may be expressed in the form 

Notice how minus two times (/x — V{x)) Q {fJ, — V{x)) in /C+,/C_ introduces the absolute value \V{x) — fj,\ into these 
differential operators. Let us denote the eigenvalues of JC± as iicj — A„, where w is a bosonic Matsubara frequency, 
and A„ is an eigenenergy of the operator —{h'^/2m)d'^/dx^ — \V{x) — //,[. The regular part of the logarithm of the 
determinant in Eq. (A12) has the sense of the free energy Fnc of an ideal gas of thermal particles 

FM - ^Indet ( ^+ ) = ^ E 1" (2 si'ili ^ 

where the regularized values of the determinants of the differential operators IC± can be obtained, for instance, by 
means of a zeta-regularization procedure [43]. Then, up to a first order in g, we find that 



—IndetG ^ = -jSFncili) + g j ^'^ j dx {Q-{x,t;x,t) + Q+{x,t;x,t)) {ipoi>o - Ptf{x; 



m)) 



= -(3Fnc{ti) + 2gj^ dr J dxpnc{x)Mo- (AU) 

Here Fn,, is the free energy of the non-ideal gas of thermal particles, and the last term in (A14) describes the interaction 
of the thermal particles with the condensate. The density of the non-condensed particles is Pnc{x) = —g±{x,T;x,T), 
and it depends, in fact, on the spatial variable only because of the translational invariance in r. For very low 
temperatures and far from the boundaries of the condensate we can use Pnc{x) — Pnc{0) because g±{x,T;x,T) is 
nearly a constant over most parts of the condensate [44] . 
We now finally obtain the effective action 

Seffitpctpo] = -PFncill) 
^ ( — S-r+K-V{x)\^bJx.T\-^ 







dr / cia;<j ■(/'o(a;,r) ^ — — A- y(a;)J Vo(a;,T) - -V'o(a;,T)V'o(a:,T)V'o(a;,T)V'o(a;,T) J- , (A15) 



where A = /i — 2gpnc{Q) is a renormalized chemical potential. We consider the Sef f (A15) as an effective one-loop 
action which involves the thermal corrections above the classical level (see Eqs. (A6)). Note that this derivation of 
the effective action does not depend on the dimensionality of the system, and is thus valid in two and three space 
dimensions as well. 

Appendix B 

We derive here the scaling exponents Eqs. (11) and (12) for infinite and small values of the coupling constant g, 
respectively, from the universal result Eq. (10) based on the Bethe Ansatz equations Eq. (4). It is possible to solve 
these equations exactly in the limiting cases g ^ and g 00. 

Let us discuss first the weak-interaction limit (/ — > 0. In order to study this free boson limit we have to rescale the 
parameters in Eq. (4): Aj = {g/Ly/'^p.j, and then take 5^0. The Bethe equations for the ground state will now 
take the form 

N 
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These equations are equations for the zeros of the Hermite polynomials. By applying the equality 



^ 1 



it becomes easy to calculate the ground-state energy of the weakly interacting gas 

In the thermodynamic limit N,L ^ oo with the density p = -j- fixed, the free energy £ of the weakly interacting Bose 

gas is 

t = hm — — = -p 

N^oo L T 

so that for small coupling constants the sound velocity is 

V = 

The substitution of this result in Eq. (10) gives the result (12). 

The case of impenetrable (hard core) bosons, also referred to as the frcc-fcrmion limit or the Girardeau- Tonks gas, 
corresponds to an infinite value of the coupling constant g — oo. The Bethe equations Eq. (4) now become 

These equations describe noninteracting particles and are easily solved, 

27r / . N+1 



and these solutions fill in the interval Ai < Aj < Ajv , where the Fermi radius Xp = Ajv corresponds to the solution 
with j = N {-Xf = Ai): 

AF = |(iV-l). 

By using the solutions of (Bl), the ground state energy of the Girardeau- Tonks gas is easily calculated to be 
In the thermodynamic limit the free energy £ of the Girardeau- Tonks gas is 



6m 

and the sound velocity at zero temperature is 

Uw 
V = — p. 

m 

Substitution of this result in Eq. (10) gives the exponent of Eq. (11). 



[1] A. Gorlitz, J.M. Vogels, A. E.Leanhardt, C. Raman, T.L. Gustavson, J.R. Abo-Shaeer, A. P. Chikkatur, S. Gupta, S. 

Inouye, T. Rosenband, and W. Ketterle, Phys. Rev. Lett. 87, 130402 (2001). 
[2] F. Schrcck, L. Khaykovich, K. L. Corwin, G. Ferrari, T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. Lett. 87, 

080403 (2001). 



18 



[3] M. Greincr, I. Bloch, O. Mandel, T. W. Hansch, and T. Esslinger, Phys. Rev. Lett., 87 160405 (2001). 

[4] K. Bongs, S. Burger, S. Dettrner, D. Hellwcg, J. Arlt, W. Ertmer, and K. Sengstock, Phys. Rev. A 63, 031602 (2001). 

[5] C.J. Pothic and H. Smith, Bose-Enstem Condensation in Dilute Gases (Cambridge University Press, Cambridge, 2002). 

[6] M.D. Girardeau, J. Math. Phys. 1, 516 (1960). 

[7] E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 

[8] C.N. Yang and CP. Yang, J. Math. Phys. (N.Y.) 10, 1115 (1969). 

[9] V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quantum Inverse Scattering Method and Correlation Functions (Cam- 
bridge University Press, 1993). 

[10] R.K. Bullough and J.T. Timonen, Quantum and Classical Integrable Models and Statistical Mechanics^ in Statistical 
Mechanics and Field Theory, Eds. V.V. Bazhanov and C.J. Burden (World Scientific Publ. Co. Pte. Ltd, Singapore, 1995) 
pp. 336-414. 

[11] D.M. Gangardt and G.V. Shlyapnikov, Phys. Rev. Lett., 90, 010401 (2003). 
[12] E.H. Lieb, R. Seiringer, and J. Yngvason, arXiv:cond-mat/0304071. 

[13] V.N. Popov, Functional Integrals in Quantum Field Theory and Statistical Physics (Reidel, Dordrecht, 1983). 
[14] T.-L. Ho and M. Ma, J. Low Temp. Phys. 115, 61 (1999). 

[15] D.S. Petrov, G.V. Shlyapnikov, and J.T.M. Walraven, Phys. Rev. Lett. 85, 3745 (2000). 

[16] N.M. Bogoliubov, R.K. Bullough, V.S. Kapitonov, C. Malyshev, and J. Timonen, Europhys. Lett. 55, 755 (2001). 

[17] S. Dettmer, D. Hellweg, P. Ryytty, J.J. Arlt, W. Ertmer, K. Sengstock, D.S. Petrov, G.V. Shlyapnikov, H. Kreutzmann, 

L. Santos, and M.Lewenstein, Phys. Rev. Lett. 87, 160406 (2001). 
[18] I. Shvarchuck, Ch. Buggle, D.S. Petrov, K. Dieckmann, M. Ziolonkowski, M. Kemmann, T.G. Tiecke, W. von Klitzing, 

G.V. Shlyapnikov, and J.T.M. Wahavcn, Phys. Rev. Lett. 89, 270404 (2002). 
[19] F. Gerbier, J.H. Thywisscn, S. Richard, M. Hugbart, P. Bouyer, and A. Aspect, arXiv:cond-mat/0211094. 
[20] S. Richard, F. Gerbier, J. H. Thywissen, M. Hugbart, P. Bouyer, and A. Aspect, arXiv:cond-mat/0303137. 
[21] K.V. Kheruntsyan, D.M. Gangardt, P.D. Drummond, and G.V. Shlyapnikov, arXivxond- mat/0212153. 
[22] E.M. Lifshitz and L.P. Pitaevskii, Statistical Physics, Part 2 (L. D. Landau and E. M. Lifshitz, Course of Theoretical 

Physics, Vol. 9) (Pergamon Press Ltd, Oxford, 1980), pp. 146-147 and following pages. 
[23] M. Naraschewski and R.J. Glauber, Phys. Rev. A 59, 4595 (1999). 

[24] U. AI Khawaja, J. O. Andersen, N. P. Proukakis, and H. T. C Stoof, Phys. Rev. A 66, 013615 (2002). 

[25] J.O. Andersen, U. Al Khawaja, and H.T.C. Stoof, Phys. Rev. Lett. 88, 070407 (2002). 

[26] D.S. Petrov, G.V. Shlyapnikov, and J.T.M. Walraven, Phys. Rev. Lett. 87, 050404 (2001). 

[27] D.L. Luxat and A. Griffin, Phys. Rev. A 67, 043603 (2003). 

[28] V. Dunjko, V. Lorent, and M. Olshanii, Phys. Rev. Lett. 86, 5413 (2001). 

[29] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[30] H.T.C. Stoof, J. Low Temp. Phys. 114, 11 (1999). 

[31] N.M. Bogoliubov, R.K. Bullough, and J. Timonen, Phys. Rev. Lett. 72, 3933 (1994). 
[32] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 
[33] S. Stringari, Phys. Rev. A 58, 2385 (1998). 
[34] V.N. Popov, J. Sov. Math. 46, 1619 (1989). 

[35] F. Dalfovo, S. Giorgirii, L.P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999). 

[36] D. Hellweg, S. Dettmer, P. Ryytty, J. J. Arlt, W. Ertmer, K. Sengstock, D. S. Petrov, G. V. Shlyapnikov, H. Kreutzmann, 

L. Santos, and M. Lcwcnstein, Appl. Phys. B 73, 781 (2001). 
[37] N.M. Bogoliubov, Zap. Nauchn. Sem. POMI 269, 125 (2000). 

[38] E.W. Hobson, The Theory of Spherical and Ellipsoidal Harmonics (Cambridge University Press, Cambridge, 1931). 
[39] V.N. Popov, JETP Letters 31, 526 (1980). 

[40] C. Itzykson and J.-B. Zuber, Quantum Field Theory (McGraw-Hill, New- York, 1980). 
[41] N.N. Bogoliubov, J. Phys. (Moscow) 11, 23 (1947). 

[42] R.K. Bullough, N.M. Bogoliubov, V.S. Kapitonov, C. Malyshev, J. Timonen, A.V. Rybin, G.G. Varzugin, and M. Lindberg, 

Theor. Math. Phys. 134, 47 (2003). 
[43] P. Ramond, Field Theory: A Modem Primer (The Benjamin Cummings Publ. Co, Massachusetts, 1981) 
[44] M. Naraschewski and D.M. Stamper-Kurn, Phys. Rev. A 58, 2423 (1998). 



